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The minimal time a system needs to evolve from an initial state to its one orthogonal state is defined as the 
quantum speed limit time, which can be used to characterize the maximal speed of evolution of a quantum 
system. This is a fundamental question of quantum physics. We investigate the generic bound on the 
minimal evolution time of the open dynamical quantum system. This quantum speed limit time is 
applicable to both mixed and pure initial states. We then apply this result to the damped Jaynes-Cummings 
model and the Ohimc-like dephasing model starting from a general time-evolution state. The bound of this 
time-dependent state at any point in time can be found. For the damped Jaynes-Cummings model, when the 
system starts from the excited state, the corresponding bound first decreases and then increases in the 
Markovian dynamics. While in the non-Markovian regime, the speed limit time shows an interesting 
periodic oscillatory behavior. For the case of Ohimc-like dephasing model, this bound would be gradually 
trapped to a fixed value. In addition, the roles of the relativistic effects on the speed limit time for the 
observer in non-inertial frames are discussed. 



Quantum mechanics acting as a fundamental law of nature imposes limit to the evolution speed of 
quantum systems. The utility of these limits is shown in different scenarios, including quantum com- 
munication', the identification of precision bounds in quantum metrology'', the formulation of compu- 
tational limits of physical systems^', as well as the development of quantum optimal control algorithms'*. The 
minimal time a system needs to evolve from an initial state to its one orthogonal state is defined as the quantum 
speed limit time (QSLT). The study of it has been focused on both closed and open quantum systems. For closed 
system with unitary evolution, a unified lower bound of QSLT is obtained by Mandelstam-Tamm (MT) type 
bound and Margolus-Levitin (ML) type bound^ '°. The extensions of the MT and ML bounds to nonorthogonal 
states and to driven systems have been investigated in Refs. 11-15. The QSLT for nonunitary evolution of open 
systems is also studied"""*. It is shown that a unified bound of QSLT including both MT and ML types for non- 
Markovian dynamics can be formulated"*. However, while this unified bound is applicable for a given driving time 
for the pure initial states, it is not feasible for mixed initial states. As we all know that decoherence and inaccurate 
operations are indispensable which may result in mixed initial states. 

Here, we shall derive a QSLT for mixed initial states by introducing relative purity as the distance measure, 
which can characterize successfully the speed of evolution starting from an arbitrary time-evolution state in the 
generic nonunitary open dynamics. Let us consider the states of a driven system in the damped Jaynes-Cummings 
model starting from a certain pure state which corresponds to a special case of our result, one may observe that the 
QSLT is equal to the driving time in the Markovian regime'". While by calculating the QSLT starting from the 
time-evolution state at any point in time which is in general a mixed state, it is interesting to find that the QSLT 
first begins to decrease from the driving time and then gradually increases to this driving time in the Markovian 
dynamics. So the speed of evolution in the whole dynamical process exhibits an acceleration first and then 
deceleration process. Additionally in the case of the non-Markovian regime, the memory effect of the envir- 
onment leads to a periodical oscillatory behavior of the QSLT. We can also focus on the widely used Ohmic-like 
reservoir spectra to investigate the QSLT for the time-dependent states of the purely dephasing dynamics process. 
We demonstrate that the QSLT will be reduced with the starting point in time for Ohmic and sub-Ohmic 
dephasing model, that is to say, the open system executes a speeded-up dynamics evolution process. While for 
the super-Ohmic environments, due to the occurrence of coherence trapping''-", we specifically point out that this 
QSLT would be gradually trapped to a fixed value, and therefore leads to a uniform evolution speed for the open 
system. We remark that the findings of those phenomena rely on our general result of QSLT for arbitrary initial 
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states. Finally, we also investigate the influence of the relativistic 
effect on the QSLT for the observer in non-inertial frames in the 
above two quantum decoherence models. 

Results 

Quantum speed limit time for mixed initial states. In the following, 
we shall consider a driven open quantum system and look for the 
minimal time that is necessary for it evolve from a mixed state to 
its final state P^+td- Under the general nonunitary quantum 
evolutions of open system, also the final state Pt_|_tj, will be 
generally mixed. One general choice of distance measure between 
two mixed states and Pi+xd fidelity F[p^,p^^^^) — 
tr ^■\/ ^PT+T^Pzy^+iD^ ■ this case of the initially mixed state p^ 
should be treated by purification in a sufficiently enlarged Hilbert 
space. And the fidelity can be written f(^ 

where the maximization is over all pure 

(Pj^)) on a larger Hilbert space that are 

^^yr-i^ defined in the smaller 
system S. However, performing the optimization over all possible 
purifications is a challenging task that will be very hard to perform 
in the general case. 

Here we consider the relative purity as a distance measure to derive 
lower bound on the QSLT for open quantum systems. The so-called 
relative purity /(t) between initial and final states of the quantum 
system is defined as^° /(t + td) = ^'■[Pi+tjjPt] /^''(p?)- To evaluate 
the QSLT, let us now characterize the derivative of the relative purity, 
f{t) = tr[p^Pi]/tr[pl). The rate of change of/(t) will serve as the 
starting point for our derivation of ML and MT type bounds on 
the minimal evolution time of an initially mixed state p-^, with the 
help of the von Neumann trace inequality and the Cauchy-Schwarz 
inequality, respectively. 

By using the von Neumann trace inequality, we begin to provide a 
derivation of ML type bound to arbitrary time-dependent nonuni- 
tary equation of the form p, = Lf (p,) . Let us consider such the evolu- 
tion. 



states 1 1//^^ 



purifications of the mixed states p^^^^^p^ 



m 



tr[pMPt)] _ tr[Lt{p,)p,] 



Hp'.) 



Hp.) 



(1) 



Then, we introduce the von Neumann trace inequality for operators 
which reads^''^^, |tr(j4iA2)| < o"i,,o"2,n where the above 

inequality holds for any complex n X n matrices Ai and A2 with 
descending singular values, o"i_i > ... > Ci „ and (T2,i — ... S a2,„- 
The singular values of an operator A are defined as the eigenvalues of 

V A'^A^^. In the case of Hermitian operator, they are given by the 
absolute value of the eigenvalues of A. We thus find 

— TTvJZ - ^ith ff, being the singular values of 

Lfipt) and qi for the initial mixed state p^. Since the singular values 
of p^ satisfy 0<p,<l, the trace norm of L^Pt) would satisfy 

\\Lt{pt)L= J2"i=i Zl"=i ^° 



En \-^n 



tr[pi) 

Integrating Eq. (2) over time from t - 
inequality 



Hp') 



(2) 



T> max 



En sr-^n 



Ttot = z + Id, we arrive at the 



\f{T + ZD)-l\Hpl)^ (3) 



where X = t^, ^ X dt.Vox unitary processes, -i'-'' equal 
to the time-averaged energy E, so the ML bound for closed systems 
can be expressed t> - |/(t + td) — l|fr(p^). 

By noting the following inequality holds ^^.„j OiQi< ffi> 



then : 



> : 



So we can simplify Eq. (3) as 



En 

[f(T + TB)-l|tr(p2 



En 



(4) 



Regarding general nonunitary open system dynamics, Eq. (4) 
expresses a ML type bound on the speed of quantum evolution valid 
for mixed initial states. 

Next we will derive a unified bound on the QSLT for the 
open systems. According to Ref 17, the rate of change of relative 
purity can be bounded with the help of the Cauchy-Schwarz 

inequality for operators, |tr(AiA2)|^ < ( AjAi ) tr( A2A2 ) . Then 



Hp.) 



tr\Li{pi)^Lt{Pf) fr(p?), since p^ is a mixed state. 



tr[p^) < 1, and we obtain 



[fit) 



tr LtiPiYltip,) 

~Hp[) 



\MpM 
Hp'.) 



(5) 



where ||i((p,) IL. = \J tr [l,(p,)'''L,(p,)] = \/Y1"= 1 

Hilbert- Schmidt norm. Integrating Eq. (5) over time leads to the 
following MT type bound for nonunitary dynamics process. 



[f(-C-hTD) 



■i\Hp'.) 



where 



. + .D 



(6) 



aj = Zjj^ V 5^ -1 ""f '^f rneans the time-aver- 

aged variance of the energy. 

By combining Eqs. (4) and (6), we obtain a unified expression for 
the QSLT of arbitrary initially mixed states in open systems as fol- 
lows. 



tqsi = max 



*[f{z + zr,)-l\Hp')- (7) 



This is one of our main results in this report. For a pure initial state 
Px=o = |'^o)(<^o|> the singular value ft = <5i,i, then 

. _ J ffj £1; = o"i < ^/y~^._^~of . We can find that the expression (7) 

recovers the unified bound of the QSLT obtained in Ref 18, which 
confirms the validity of our results. So Zqsi formulated in (7) defines 
the minimal evolution time for arbitrary initial states. 

The speed of evolution in the exactly solvable open system dyna- 
mics. In order to clear which bound on the speed limit time Zqsl can 



be attained and tight, we must compare aiQj and y _ j 

For instance, when ""ift < ^/'y~^._,"of , the ML type bound 

provides the tighter bound on the QSLT. Next we shall illustrate the 
application of the QSLT to the quantum evolution speed of a qubit 
system in two decoherence channels. A generally mixed state po of a 
qubit can be written in terms of Pauli matrices, whose coefficients 

define the so-called Bloch vector P(,= - {l + Vj^ax + Vyay + v^a^), 
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where / is the identity operator of the qubit, (Tj- (fc = x,y, z) is the PauU 
operator, and v= (v;c,v^,Vz)e_B : = {vejf'; || v|| < l}. 

We firstly consider the exactly solvable damped Jaynes- 
Cummings model for a two-level system resonantly coupled to a 
leaky single mode cavity. The environment is supposed to be initially 
in a vacuum state. The nonunitary generator of the reduced 

dynamics of the system is L((p() =y, ^cr — P(t7+ — -0"+ t7_P( — 

-p,o"+t7_), where a± = ± iOy are the Pauli operators and yt is 

the time-dependent decay rate. In the case of only one excitation in 
the whole qubit-cavity system, the environment can be described 
by an effective Lorentzian spectral density of the form 
1 y ). 

J(co) = = where a is the width of the distribution, 

' 2n{ojo-mf + ?.^ 

cUq denotes the frequency of the two-level system, and y'o is the coup- 
ling strength. Typically, weak-coupling regime is, ). > 2yo. where the 
behavior of the qubit-cavity system is Markovian and irreversible 
decay occurs. The strong-coupling regime is, a < ly^, where 
non-Markovian dynamics occurs accompanied by an oscillatory 
reversible decay. The time-dependent decay rate is then explicitly 



given by = 



ly^X sinh{dt/2) 



. , , , . , , , , , . with d = \ )? - 

d cosh( A/2) + ). sinh(dt/2) V 

The reduced density operator of the system at time t reads 



-2yoA. 



2-{l-v,)p, 



(l-^z)p, 



where pt = e Jo * "' . 

For the generally mixed state of a qubit, QiOi + p20"2 is always 



less than y erf -|- o'2> so we reach the result that the ML type bound on 
the QSLT is tight for the open system. The unified expression (7) 
proposed for the mfaced initial states can demonstrate the speed of the 
dynamics evolution from an arbitrary time-dependent mixed state p^ 
to another P^ + ^d ^7 ^ driving time Tq. We examine the whole 
dynamics process where the system starts from the excited state, 

= — 1 and = Vy = 0. Figs. 1(a) and 1(b) show the QSLT for a 
time-dependent mixed state as a function of t in the Markovian 
and non-Markovian dynamics process, respectively, where = 1. 
The QSLT can initially decrease to a minimum, and gradually 
reaches to the driving time t^, in the Markovian regime. While for 
the non-Markovian regime, the QSLT decreases to a minimum in the 
beginning of the evolution, then occurs a periodical oscillatory of the 
time T. That is to say, in the Markovian regime, the evolution of the 
qubit first exhibits a speeded-up process for t < and then shows 
gradual deceleration process for t > t^. In contrast, the speed of 
evolution for the qubit in the non-Markovian dynamics process 
complies with an interesting but reasonable periodical oscillatory 
behavior. This is a newly noticed phenomenon. 

The above behavior can be explained by evaluating the QSLT for 
the qubit to evolve from p^ to p^+ij,. 

!(pT-pz + rJ(l-2p.)| 



dt 



(9) 



For the Markovian regime y, = yo, the value of ]pj ] can be given by 
yje"'"', then the speed limit time is simplified as 
tqsi = td| 1 — 2e^''°^|. So the appearance seen in Fig. 1(a) depends 
only on the decay of the excited population pj = e^' "^ for the time- 
dependent state p^, and the critical time = — In 2. In comparison, 

"/o 

the oscillatory behavior shown by Fig. 1 ( fo) in the the non-Markovian 
regime appears as a consequence of the oscillatory time dependence 
of the decay rate )',. 




Figure 1 | The QSLT for tlie damped Jaynes-Cummings model as a 
function of the initial time parameter r. Here we choose the excited state 
as the initial state, Vz = = Vy = 0. (a) the Markovian regime, 

}'o = 0.1a. (b) the non-Markovian regime, yo = 10/. The red (dashed) line 
represents the decay of the excited population p^. Parameters are chosen as, 
coo = 1> / = 1 and = 1. 

As an additional application, when the initial state takes the form 
of general mixed state, the speed of the qubit dynamics process in the 
damped Jaynes-Cummings model can also be studied, the results are 
presented by the dashed red line and dotted blue line in Fig. 2. We can 
find that in the Markovian regime, the QSLT initially reduces to a 
minimum and then increases to a maximum. Next, the QSLT would 
gradually decrease in the following dynamics process. This can be 
understood that the evolution of the whole dynamics process starting 
from a mixed initial state exhibits an eventual acceleration process 
after a speed-up and deceleration process in the Markovian regime. 

In what follows, we consider a spin-boson-type HamUtonian that 
describes a pure dephasing type of interaction between a qubit and a 
bosonic environment. It is worth stressing that this qubit-environment 
model admits an exact solution^''^''. There exists no correlations 



(0 
O 




Figure 2 \ The QSLT for the general mixed states. In the Methods section, 
we have obtained the QSLT for the qubit to evolve from to P^+td ''^ 
damped Jaynes-Cummings model with the specific form Tq"^. From the 



it is easy to show that Tqj^ is relate to and the 
vl + vi of the initial state. So this figure describe the 



expression of Tq°^, 
coherence C(po) = 
speed of the qubit dynamics process starting from the general mixed initial 
states in the damped Jaynes-Cummings model in the Markovian regime 
(yo = O.l/l). Parameters are aiti = \, X = \ and To = 1. 
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between the system and the environment at t = 0. Furthermore, the 

environment is initially in its vacuum state at zero temperature. The 

nonunitary generator of the reduced dynamics of the system is Lt(P() 

= JtiOzPt^z ^ Pt)/^- By considering that the bosonic environment 

operator is simply a sum of linear couplings to the coordinates of a 

continuum of harmonic oscillators described by a spectral function 

, , , / Hffl \ 1 — cos cot 
dw]{(j}) coth „ I . Here, we 



then 



2knT 



suppose that the spectral density of the environmental modes is 



Ohmic-like /(cu) = - 



, with oj^. being the cutoff frequency 



and 71 a dimensionless coupling constant. By changing the s-para- 
meter, one goes from sub-Ohmic reservoirs (s < 1) to Ohmic (s = 
1) and super-Ohmic (s > 1) reservoirs, respectively. For zero temper- 
ature, t > 0 and s > 0, the dephasing rate can be obtained 

cos[(s — 1) arctan(a)cO]r(5— 1) 



1- 



where r(s — 1) is the 



(1+C0?f2) 

Euler Gamma function. Taking the limit 5^1 carefully, one also finds 
y{t,s= l) = r]\n(^l+o)lt^). The time evolution of the reduced density 
matrix of the qubit satisfies 



1 



l+Vz 
{Vx + iVy)qt 



(1-v. 



(10) 



where qt — e^" . 

In this Ohmic-like dephasing model, the QSLT of a qubit can also 
be described by ML type bound. In the dephasing evolution, by 
considering an arbitrary mixed state evolves to another Pj+^j, 
under a driving time Tq, the QSLT can be calculated 



C(Po)^^^|<?T?T+ic-q?| 



1 

td J 



(11) 



\qi\dt 



With this, it is easy to show that Tqsl is independent of v^, and not 
only it relates to the dephasing rate of the Ohmic-like environment 
but also to the coherence of the initial state po under a given driving 
time Xr>- Fig. 3 presents the results of our analysis for Tqsl in the 
Ohmic-like dephasing process with different C(po)- We observe that, 
for the same driving time = I, the lager coherence of the initial 
state can decrease the speed of evolution of a quantum system, and 
thus demand the longer QSLT. By choosing the 5-parameters which 
correspond to Markovian regime^', the speed limit time can be 
rewritten as Tq^^ = T£)C(Po)"'''Jt- Hence, for a given initial state po, 
the speed of evolution in the dynamical process is determined by the 
decay rate of the coherence for the mixed state p^. Due to the 
specific form of the spectral density for Ohmic-like dephasing 
model", in the case of zero temperature, the qubit dephasing q.^ will 
predict vanishing coherences in the long time limit for s £ 1. On the 
other hand, for s > 1 the qubit dephasing q^ wiU stop after a finite 
time, therefore leading to coherence trapping, as shown by the inset 
of Fig. 3. This is important for quantum information processing since 
coherence of the quantum state is preserved. Other observations in 
Fig. 3 are as follows: The open system executes a speeded-up dynam- 
ical evolution process in the Ohmic and sub-Ohmic dephasing mod- 
els. But for the super-Ohmic dephasing model the qubit firstly 
exhibits a speeded-up dynamical process before a finite time, and 
then complies with an uniform evolution speed after this finite time. 

Relativistic effect on quantum speed limit time. If the observer for a 
quantum system in a uniformly accelerated frame with acceleration 
a, the relativistic effect should be taken into account^" Here we 
briefly present the study of the relativistic effect on the QSLT. Owing 
to the relativistic effect, the coherence of the changing initial state 

v^ + Vyj, and becomes much less than that of C(po). 



turns mto cos r 



-C(p„)=1,s=2 
C(p„)=0.5, s=2 
• C(p„)=0.25, s=2 
-C(pJ=1,s=1 




Figure 3 | The QSLT for the Ohmic-like dephasing model as a function of 
the initial time parameter r with different coherence C{/>o)- The inset 
shows the decay rate of the coherence for the mixed state Pi with C(po) = 
L Parameters are chosen as, co^ = 1, f; = 1 and = L 

So the relativistic effect can increase the speed of evolution of a 
quantum system in the purely Ohmic-like dephasing channels. 
The parameter r above is defined by cosr = (g-^'"'''^''" -|- 1)""^, c the 
speed of light in the vacuum, and tu the central frequency of the 
fermion wave packet. But for the damped Jaynes-Cummings 
model, in spite of the weaker coherence of the initial state brought 
by the relativistic effect, the larger excited population 
(l-l-Vz)cos^r 

1 ^ m the changing initial state can also be acquired. 

As well as the QSLT mainly depends on the population of excited 
state under a given driving time in the amplitude-damping 
channels''\ the relativistic effect would slow down the quantum 
evolution of the qubit in the damped Jaynes-Cummings model, 
therefore leads to a smaller QSLT. 

Discussion 

We have derived a QSLT for arbitrary initial states to characterize the 
speed of evolution for open systems. Some novel phenomena are 
observed. For the damped Jaynes-Cummings model, we have 
obtained that the speed of evolution, where the qubit starts from a 
pure initial state, exhibits an acceleration first and then deceleration 
process in the Markovian regime, and shows a peculiar periodical 
oscUlatory behavior in non- Markovian regime. Moreover, in the case 
of the purely dephasing environments, the QSLT would be gradually 
reduce to a fixed value for the super-Ohmic dephasing model, and 
hence leads to a uniform evolution speed for the open system. Our 
results may be of both theoretical and experimental interests in 
exploring the speed of quantum computation and quantum informa- 
tion processing in the presence of noise. 

Methods 

Quantum speed limit time in open system dynamics. Let us consider the exactly 
solvable damped Jaynes-Cummings model and the Ohimc-like dephasing model. By 
considering the initial mixed state 



its time evolution reads 



1 + v, 



1- 



(12) 



(13) 



Then we will demonstrate the QSLT for the open system dynamics process from an 
arbitrary time-dependent mixed state p-^to another Pr+io by a driving time t^- Due to 
the unified expression for the QSLT of arbitrary initially mixed states in our work, we 
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first calculate the singular values of and Ldpt), respectively. For p^^, the singular 
values — 1, 2) are 



(14) 



02 = -^[^ + , 



where v'x — Vxy/p^, v'y — Vy-Jp^ and — 1 — {)- — Vz)px in the damped Jaynes- 
Cummings model. While for the Ohimc-like dephasing model, v'x — Vxqz, v'y — Vyq^ 
and v'z — Vz- And the singular values o"; of I[{pf) can be expressed as 



4(l-vj' + - 



(15) 



So we can obtain that q\G\ + Qi^^i is always less than \J <^\ + ^n*^ reach the result 
that the ML type bound on the QSLT is tight for the open system. The QSLT for the 
qubit to evolve from p^ to P^+t;^ in the damped Jaynes-Cummings model can be 

obta ned t?£ - , , ' , 

QSL 1 + / 1 \ Jj. 



.where [/(i + Td) - 1| tr(pj) = 



\tr{p,+,,P,)-tr{p',)\ = l\{l-v,){p,-p,+,„)ll-(l-v,)p,]+(vl + vf) 

[y/pzpT + To ~pT ) I - With this expression, it is easy to show that Xqsl is not only relate to 
but also to the coherence of the initial state C(po)- With regard to the QSLT for the 
qubit to evolve from p^ to p^_|_^^ in the Ohimc-like dephasing model, we can obtain as 



rT + Ti, 



Quantum speed limit time in non-inertial frames. If several detectors or observers 
are located in different inertial or non-inertial frames, or in curved space-time, the 
relativistic effect should be taken into account. So in this subsection, we may 
investigate the influence of the relativistic effect on the QSLT in the damped Jaynes- 
Cummings model and the Ohimc-like dephasing model. To start with, let us consider 
the generally mixed state po in Eq. (12) of a qubit held by Bob in a non-inertial frame. 
Now we assume that Bob moves with a uniform acceleration and takes with him a 
detector for the qubit. From the perspective of Bob, under the single-mode 
approximation, the Minkowski vacuum state |0)jvi and the Minkowski only excited 
state 1 1 )m, in a uniformly accelerated frame with acceleration a are expressed as^°, | 0)m 
= cosrjO)/ 10)77 + sinr|l)7 |l>//, and \1)m ~ \ The parameter ris defined by cosr 

= I" + l)-"\cthe speed oflightin the vacuum, and tjj the central frequency of 

the fermion wave packet. And the subscripts I and II represent the states related to the 
Rindler- region -I and -region-II respectively. Here r ranges from 0 to nIA for 0< a < 

Due to the presence of a Rindler horizon, Bob is forced to trace over a causally 
disconnected region-II of space-time. Accordingly, the density matrix of the initial 
state Po becomes 



1 / (l-Fv^)cos^ 



2 \ (vx -\- iVy^ cosr 2sin^r-|- (1 — V2)cos^r 



(16) 



Then the time evolution of the reduced density matrix pjf in the damped Jaynes- 
Cummings model and the Ohimc-like dephasmg model can be written as follows, 
respectively. 



Pn 



2 — [2sin^r+ (1 — Vz)cos^r]p, (v, — iv^))cosr,/p7 
(vj + iVy)cosr^/pi [2sin^r+ (1 — V2)cos^r]p, 



and 



Pn 



1 / (l + v^jcos^r [vx — iVy))cosrqi 

2 \ (vx-\-iVy)cosrqi 2sin^r-|- (1 — V2)cos-^r 



(17) 



(18) 



According to the relativistic effect, the QSLT for the qubit to evolve from po to p,^ in 
the damped Jaynes-Cummings model can be shown as the form t^^^ — 

iJo (giai + eioijdt ^ ' 41 \ '7 

( 1 + ) cos^ r + [2sin^ r + ( 1 - Vz ) cos^ r] ^ (pi - 1 ) - ( 1 + ) cos V [2sin^ r + 



(l — Vz)cos' r]p,\, and (cos2r — 3 + 2cos^rVz)'^ + cos^r|^vJ + j ^ pt = 

gi(7i + f)2C2- while for the case of the Ohimc-like dephasing model, the QSLT can be 
cosr + vj 1 — ijt 
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